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fH , Abstract 

^ . Recent data on inclusive B decays into charmonium states taken by the CLEO 

^ ■ collaboration suggest that a sizeable fraction of the nonleptonic decays may consist 
of the J/ipAp three body final state, corresponding to a distinct enhancement 

^ I in the inclusive J/V' momentum distribution. The kinematical boundary of this 



structure corresponds to the case where the J/il^ recoils nearly monoenergetically 
in the B rest system against a partner having mass of ~ 2 GeV. This may allow 
the observation of a A — p bound state near or just below threshold; i.e., strange 
baryonium, even if the production rate is small. Using a phase space approach 
to the B meson decay, we study the J/ip momentum distribution and the effect 
of baryonium formation. We also discuss the possible observation of pentaquarks 
and hadronically-bound J/ip by the observation of monoenergetic baryons in these 
decays. 
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1 Introduction 



One of the most important new elements that quantum chromodynamics brings to the 
strong interactions is the potential existence of new bound states and resonances be- 
yond the conventional mesons, baryons and nuclei. These include exotic states such as 
gluonium {gg, ggg), hybrids {qqg) and dibaryons such as the H (udsuds), etc. 

A particularly sensitive measure of QCD binding is nuclear-bound quarkonium, i.e., 
a (c — c) or (b — b) pair bound to a nucleus. In contrast to the usual strongly interacting 
few-body systems, the hadrons in the (cc) (uud) meson-nucleon system have no valence 
quarks in common. However multiple-gluon exchange (the QCD van der Waals force) 
may lead to the formation of bound states [|1| . Using the fact that the characteristic size 
of the -J/ip is much smaller than the typical hadronic scale {^qcd), a QCD calculation 
based on a multipole expansion was performed in [0, leading to nuclear-bound quarko- 
nium binding energies of order 8 — 11 MeV. In addition, attractive interactions may 
also arise from the J/ipp D Ac interaction ||^. Long-range exchange interactions 
may arise in second order perturbation theory if there are significant nonvalence [QQqq) 
Fock states in the quarkonium wavefunction or intrinsic charm [uudcc] Fock states in 
the nucleon. There is thus a possibility of a strong J / ip-micXeow resonance or bound 
state. 

From the experimental point of view, there is evidence for the state = \qqqss> 
(where q are u or d quarks), reported by the SPHINX collaboration [Q, and there are 
proposals for a search for the = \qqqcc> with the same setup and at Fermilab 
as well [0. 

Another interesting type of exotic state that could arise in QCD are bound states 
of baryons and antibaryons (baryonia) such as p — p and A — p, and pentaquarks, Pq, 
baryons containing five quarks like uudcs or uudcs Baryonia may have binding 

energies as small as a few MeV as is characteristic of nuclear interactions, but they are 
unstable and decay into mesons via q — q annihilation with possibly narrow widths. Such 
states have been subject of intensive searchf^]. In particular, experiments have looked for 
strange baryonium in hadronic collisions, typically in Kp ^ App reactions at laboratory 
energies of 8.25 GeV [§, 18.5 GeV [10| and 50 GeV |11|. The results concerning the 
width of these states have been somewhat controversial. Some experimental searches 



r|] found r ~ 20 MeV whereas others [0, ^ found F ^ 150 MeV. The existence of 
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narrow structures has never been established or completely ruled out. 

Searches for the pentaquarks are in progress at Fermilab. Data from the collaboration 
E-791 are currently being analyzed, and another experiment, the E-781, is accumulating 
events with charmed baryons, which may determine the existence of certain pentaquarks 

In this paper we will show how to look for QCD exotic states in a B-factory. In 
particular, we will show how the J/ipAp decay channel can be used to search 

for the {J/ip—p), [J ftp — A) and (A — p) bound states. Also the appearance of two- 
body structures in the momentum distribution of the Ac near p = 0.47 GeV/c in the 
B^ —>■ Ac + X decay channel may provide a signature for the nude's pentaquark. 

2 B~ J/ijjKp decay mode 

Among the numerous decay modes of the B mesons, the B — > charmonium + X channels 
have drawn particular attention. Recent measurements of these decays performed by the 
CLEO collaboration |T^ have found the following inclusive branching fractions: (1.12 ± 
0.04±0.06)% for B J/^X, (0.34±0.04±0.03)% for B ^ i;' X and (0.40±0.06±0.04)% 
for B XciX. Crossing this information with the measured |jl5| branching fraction 
(2.5 ± 0.4)% for B ApX, one can infer that a significant fraction of the charmonium 
decays correspond to the three-body decay: 

B-^J/iPAp (1) 

This process, ideal for our purposes, is depicted in Fig. la. In first place, the e^e+ BB 
environment is much cleaner than in hadronic collisions, and the four momentum of the 
B is precisely determined. Moreover these decay products are very massive and only a 
small amount (~ 128 MeV) of the initial energy is converted into kinetic energy. The 
decay particles move slowly in the B rest frame, interact strongly and can thus form 
bound states just below threshold or resonances just above threshold. 

The CLEO collaboration has already provided a measurement of the inclusive mo- 
mentum distributions of the J/ijj and ip' in B decays [T^]. If the A — p is formed as a 
bound state, the J/ip will be produced with a momentum just above 0.56 GeV/c (in the 
B rest frame). Similarly, if the J/ip — A is bound, the p will be produced at a momentum 
above 0.45 GeV/c. Finally, if the J/tp — p is formed, the A will have momentum slightly 
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larger than 0.48 GeV/c. If the state is below the J /if) — p mass threshold it may 
be observed in the decay as a narrow peak in the A momentum distribution at 0.48 
GeV/c. Another interesting decay channel is shown in Fig. lb, where the proton and 
the meson may be bound forming the nude's pentaquark. The conjugate process 
is shown in Fig. Ic. In Fig. Id we illustrate the B^ Dj pp decay. Here the proton 
and the D~ meson may be bound forming the uudcs pentaquark. This state, although 
CKM suppressed, is interesting because it is more strongly bound than uudcs [7|. 
It thus seems possible, with high enough statistics, to identify two-body structures in 
these spectra, direct signatures for the formation of exotic QCD bound states. 



In Fig. 2 we reproduce the B J/ipX spectrum presented in |Tj], in which the 
feed-down modes B ijj' X and B XaX have been subtracted. This procedure 
eliminates the J/ip^s coming from decays of excited charmonium states. The number of 
remaining points is not large, but we can nevertheless clearly see a distinct secondary 
bump in the lower momentum region corresponding to the B J/ippA decay. 

A baryonium state such as A — p has the same quantum numbers as an excited kaon. 
We can distinguish a nuclear-type A —p bound state by the fact that its mass should 
be within lO's of MeV of the threshold. Thus B decays offer not only the possibility of 
observing baryonium resonances ||, ^, |ll|, |12| , but also bound-state baryonia such as 
the A— p with mass below the threshold ( 2.05 GeV). The (A— p) state recoils against the 
J/ip with momentum (in the B~ c.m.s.) slightly beyond the upper limit (0.56 GeV/c) 
of the phase space for the three body decay B ^ J/tp Kp. 

If there is a narrow peak at pjj^ ~ 0.56 GeV/c, the mass Mx of the J/ip decay 
partner can be immediately inferred from the simple two body phase space formula: 

with 

A (x, z) = x'^ + y'^ + z"^ — 2xy — 2xz — 2yz (3) 

In the above expression Mb is the B meson mass (= 5.279 GeV) and mj/^ = 3.097 GeV. 
This simple equation can be solved for Mx giving Mx — 2.05 GeV, which is exactly 
equal to the sum of the A and p masses. A peak above 0.56 GeV would be an indication 
of a Ap bound-state formation. If the binding energy of the A — p system is Amxp — 10 



4 



MeV then 



^P%, = ^<f^ (4) 

giving Apj/^ = 21 MeV/c beyond the upper hmit of the allowed kinematical region for 
B — s> J/ipAp. Given that the estimated experimental resolution is of Ap =1.7 MeV 
|T6| , it seems possible to detect such peak. 

In Ref. the ijj' spectrum is also presented. The statistics are still small, and 
there is no evidence of any structure. This is in agreement with our phase-space based 
considerations: because of its larger mass (~ 3.686 GeV), the ijj' decay partners can have 
masses of at most 1.593 GeV. This is not enough to produce two baryons and therefore 
a near-threshold decay is excluded. 



3 Modified phase-space analysis of 5-decay 

Quantitative predictions for exclusive decay amplitudes such a.s B —>■ J/ tpAp or strange 
baryonium formations B J/iplAp] involve all of the complexities and uncertainties 
of QCD hadronization at the amplitude level and at low relative momentum where the 
gauge interactions are strongest. Nonrelativistic QCD (NRQCD) has brought significant 
progress to this field. The underlying weak decay b ccs provides a guide to the flow 
of the heavy quarks. Very recently [|T^ the J/ip momentum distribution in the decay 
B ^ J/ip + X was calculated for the first time and compared to data. It is interesting 
to observe that the obtained spectra agree with data, but fail in the low momentum 
region, exactly where we may expect some enhancement due to the exclusive channel 
considered here. 

In general, hadronization will distort the "bare" quark momentum distribution given 
by the weak matrix element. Another source of uncertainty is the momentum distri- 
bution of the initial quarks inside the B meson, which was shown in [1^ to strongly 
affect the spectrum of final particles. In order to keep our approach simple, we shall 
assume that the J/ip spectrum measured over the whole available phase space gives a 
good indication of how the complicated underlying dynamics modifies the pure phase 
space picture ( = const) , thus providing a good representation of the squared 

matrix element. 
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In what follows we describe our modified phase space analysis of B decay. Following 
the standard description of the kinematics JTB] the decay rate, at the hadronic level, is 
given by: 



E 1-^ r (2^)'-' n iP- - P-'^ - pS) 6{V~s~E,-E,- E,) (.5) 

In the above expression the indexes 1,2 and 3 refer to the J/ip, A and p respectively. Ei 
is the energy of the i^'^ particle and ^/s is the center of mass energy. We take |A^p as 

— ^ const -/(pi) (6) 

where the function f is given by 

fix) = X ■ exp l^-^-^-p^y {1 - x) (7) 

where x = pi / Punix, xq = 2.0 / Pmax, o"o = 0.92 / pmax and pmax is the maximum momen- 
tum allowed for the J/ip. In our case pmax = 1-95 GeV and the function f reproduces the 
data presented in [Q. The normalization of / is absorbed in the constant in eq. (^). 

In the rest frame of the decaying particle (p = and ^/s = Mb) the integration 
over the variables 2 and 3 can be then easily done and we obtain the J/ip momentum 
spectrum: 

f- = ^m)-X'/'{Mls„ml)^'^'{s2,ml,ml) (8) 
dpi El S2 ^ ^ ^ 

where we have omitted all constant factors (since the normalization is free) and S2 = 
+ mf — 2MbEi. Inserting in eq. (P) the values of the masses, we obtain a spectrum 
which is adjusted to data and plotted in Fig. 2 as a solid line. As can be seen, the 
position of the maximum given by the simple modified phase space model coincides 
with the maximum of the bump. In Fig. 2 we also show (with a dotted line) the 
fit of data with eq.(^. In the region pj/^ < 0.56 GeV/c this curve is an estimate 
of the background, i.e., all processes other than eq. (|l|) in which a J/ip is produced. 
Subtracting the integral of the dotted curve (over < pj/^ < 0.56 GeV/c ) from the 
integral of the solid line we can obtain an estimate for the branching fraction for eq. (0) 
~ 0.04%. 

For a baryonium binding energy of Am\p = 10 MeV, the exact position of the peak 
in the J/ip distribution is given by eq. (0) with Mx = rn^ + nip — Am\p. We assume 
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that the width of this state is equal to its binding energy. A final condition must be 
satisfied by this bound state: it will be observed only if the strength of the peak is much 
larger than the background. For illustration we have estimated this strength to be 10% 
of the whole distribution. This produces a clear signal, which is shown in Fig. 2 with a 
dashed line. More details are given in the next section. 



4 The strength of the peak 

In this section we explain how to estimate the strength of the quasi-two-body decays cor- 
responding to the area below the (dashed line) peak in the J/ ip momentum distribution. 
The process 

B- — ^J/ijAp (9) 

can be understood, at the quark level, as the reaction 

b +u — > ccud suud (10) 

characterized by a cross section, a, which cannot be rigorously calculated. Since we are 
always treating the case where we have a J/ip , we can coalesce the two charmed quarks 
into one particle and say that the effective reaction is 

b +u — > J/i/j + udsuud (11) 

where the six quarks system has invariant mass mg. Following the spirit of the ap- 
proximations made in the last section, the differential cross section a can be written 

as 



da = const ■ dRj (12) 

where the constant factor represents the matrix element squared and R7 is the phase 
space factor for the seven final particles {J/ip udsuud), which, with the help of a 
recursion relation ||18|, can be written as a function of mg : 

dRj = ^^R,{m,) (13) 
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In the center of mass system, from energy momentum conservation we have 



Mb = + + V^j/^ + ^i (14) 

Solving this relation for pj/^, differentiating with respect to mg, substituting in eq. ( pTSD 
and inserting eq. (|13|) into eq. ( ]I2| ) we obtain 



, / \ 3/2 

"0" nz / "^6 



where a nonrelativistic approximation is made and some immaterial constants are dropped. 
With the additional assumption that the six body phase space is approximately one- 
dimensional and taking constituent masses for the quarks (~ 300 MeV) we can evaluate 
i?6 (almost) analytically. As expected the curve Eq. (|1^) has a bell shape. Rq starts 
from zero when m, is just the sum of the quark masses and then increases monotonically. 
The multiplying factor in Eq. (|15D is a decreasing function of m, and goes to zero when 
it is equal to Ms — mj/^. This cross section can be integrated over different mg intervals. 
The bound-state region corresponds to m, smaller than the A — p threshold, i.e., 

mp + m\ — e < m, < mp + m\ (16) 

where e is the binding energy. The continuum region is given by 

rrip + mA < me < Mb - mj/^ (17) 

The region 

me < mp + m^ — e (18) 

refers to processes in which the B mesons decays into pions and a kaon with too small 
invariant mass. Such events do not come from the baryonium decay and are thus unin- 
teresting. This part of the phase space will be ignored. 

The probability of forming a bound state is then given by 



_ Jmp+mA-e drrif, qn 



We use e = 10 MeV. Performing the calculations we find that p = 10%. 
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On the other hand, looking at Fig. 2, we observe that the probabihty of forming a 
bound state may be identified with the ratio between the area below the (dashed line) 
peak and the total area (solid plus dashed lines), i.e., the differential branching fraction 
integrated over pj/^ : 

where ^pj/ip given by Eq. (^). Knowing p from eq. (|l^), we replace it in the above 
expression and determine the right hand side, i.e., the strength of the peak, which, with 
the parameters stated above, is 10%. 

The same procedure can be applied to the other interesting bound states in the decay 
mode (|I|) as mentioned at the end of the introduction. The {J/ip — p) and (J/V" ~ A) 
may be observed as peaks in the momentum distribution of the A and p respectively. 
These spectra are computed with the help of eq. (|). Now, pi refers to p or to A and 
/(Pi) = I5 since we do not know how to modify the phase space in these cases where 
no data are available. The strengths of the bound state peaks are computed in the 
way described above (with e = 10 MeV and Anij/^p = Amj/^\ = 10 MeV), and the 
results are p = 10.7% and p = 11% for the A and p spectra respectively. These spectra 
(with corresponding peaks) are shown in Fig. 3 in arbitrary units. 



5 Conclusions 

The high luminosity of the CESR storage ring and the 5-factories now under construc- 
tion will provide a high statistic sample of B decays. As we have discussed in this 
paper, the decay of the B may lead to the production of exotic QCD states such as 
quarkonium-baryonium bound states, pentaquarks, and strange baryonium. The nearly 
monoenergetic production of the recoil hadron in the rest system of the B could provide 
a sensitive and clean signal for these states. 
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Figure Captions 



Fig. 1 Exclusive hadronic B decays which may lead to the formation of bound states, 
la) hadronically bound J/ip or strange baryonium; lb) nude's pentaquark; Ic) 
uudcs pentaquark; Id) uudcs pentaquark. 



Fig. 2 The J/i/j momentum spectrum measured by the CLEO Collaboration |T^ after 



the subtraction of the feed-down modes. The solid line is the result of a modified 
phase space calculation for the B Jip Ap three-body decay, eq.(|). The dotted 
line shows the fit of data, eq.(^, used to estimate background contributions. The 
dashed line represents the expected peak if there is a narrow strange baryonium 
with binding energy of 10 MeV and width 10 MeV. 

Fig. 3 The momentum distribution of the p (solid line) and A (dashed line) produced 
in the B decay eq.(p. The broad curves show the three-body decay spectra of 
the p and A calculated with eq.(P). Just beyond the kinematical limits of these 
two curves we show the peaks which would be signatures of the below-threshold 
(J/ip — A) (sohd line) and (J/i^ ~p) (dashed line) bound states. Their binding 
energies and widths are assumed to be e = 10 MeV. These numbers determine the 
position and width of the bound state peaks. Their strengths are given by eqs. 
dT^) and (gup. 
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